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We examine the possible resolution of the recently observed polarization anomaly in B°(B°) — > 
(f)K*° (K*°)-decay within R-parity violating ($p ) SUSY. We show that a combination of the su- 
perpotential trilinear $ p -interactions, with the couplings A', and the soft SUSY breaking bilinear 
]fi p sneutrino-Higgs mixing, proportional to p 2 , can potentially generate the effective operators 
with the chirality structure necessary to account for this anomaly. However, we demonstrate that 
the existing experimental data on B s — » fi + fi~ -decay lead to stringent upper limits on the Wilson 
coefficients of these operators, which are about two orders of magnitude below the values required 
for the resolution of the B-decay polarization anomaly, and, therefore, it can hardly be explained 
within the Ijtp SUSY framework. 

As a byproduct result of our analysis we derive new limits on the products of the soft bilinear 
and the superpotential trilinear ^-parameters of the form jl 2 \' . 
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I. INTRODUCTION 



Now it is widely recognized that B-mesons offer powerful means for testing the standard model (SM) and probe 
physics beyond its framework. Recently, remarkable progress has been achieved in experimental and theoretical 
studies of B-physics. One of the most important experimental results of the last years in this field was, certainly, the 
discovery of CP violation in the B-system. The running B-experiments [J at BABAR, BELLE, CDF, DO and CLEO 
have also collected a large statistics on various decay modes of B-mesons some of which seem to be quite challenging 
for the SM. 

The BABAR [2] and BELLE [3] Collaborations reported experimental data on B-meson decay to a pair of light 
vector mesons: B — ► VV where V — p,4>,K*. An intriguingly large transverse polarization fraction comparable to 
the longitudinal one has been observed in the B°(B Q ) — > <pK *°(K *°)-decay channel. This result has been recently 
confirmed by the CDF collaboration Q as well. This polarization anomaly is hard to be explained within the SM 
and may indicate some new physics. As is known, the SM predicts for the helicity amplitudes of B° — > (j)K*° the 
following ratios [EH| : Hoo ■ H— : H ++ ~ 0(1) : 0(l/m,b) : 0(l/m%), where H 00 corresponds to the final vector 

mesons in the longitudinal polarization state, while H + ±,H in the transverse positive and negative one. This SM 

result is in an obvious disagreement with the BABAR 2], BELLE Q and CDF Q observations, demonstrating that 
|ff ++ ±H__| a »|ffoof. 

In the literature various efforts have been undertaken to account for the polarization anomaly from the view point 
of the SM 0, 0] and in various scenarios of new physics beyond the SM [!,[§]. In Ref. [§] a model independent analysis 
of the B-decay polarization anomaly has been carried out on the basis of the effective Lagrangian approach. Two sets 
of effective AB = 1 operators necessary for the resolution of this anomaly have been identified. In addition from the 
experimental data 2, 3] the corresponding values of their Wilson coefficients have been determined. These two sets of 
operators have the following chirality structure: (i) (1 — 75) ® (1 — 75), cr(l — 75) ® a(l — 75) and (ii) (1 + 75) <E> (1 + 75), 
(7(1+75)^(7(1 + 75). 

In the present paper we use this model independent result in order to examine the ability of R-parity violating 
SUSY ($ p SUSY) to resolve the above discussed polarization anomaly in B a (B°) — > cf)K* (K* a )-deca,y. In Sec. 2 we 
specify the effective AB = 1 operators satisfying the polarization Anomaly Resolution Criteria (pARC). In Sec. 3 
we determine these operators in the context of Ifi p SUSY and derive their Wilson coefficients. In Sec. 4 we study 
experimental limits on these Wilson coefficients from the existing B s — > fi + fj," data and discuss the compatibility of 
the pARC with these limits. 



II. CRITERIA FOR RESOLUTION OF THE POLARIZATION ANOMALY 

The effective Hamiltonian Ti describing B° — ► <pK*° with AB = 1 can be written in the form: 

14 44 

n AB =i = ni K Li + <l=i = -§ E «M ■ Oi{n) + -§ £ c^) • o^) + h. c . , (i) 

where Ci(/x) are the Wilson coefficients evaluated at the renormalization scale /i ~ mj,. The first 14 terms correspond 
to the penguin-dominated SM contributions T~Cab=i listed in Ref. [lOj, the last 30 terms Wab=i appear in the presence 
of new physics (NP). In Ref. [9( it was shown that out of the 30 NP-operators only the following operator set 

O15 = S a P R b a ■ S(jP R Sf3, Ol6 = S a P R bfj ■ SpPRSa, (2) 

On = s a P L b a ■ spPLSp, Ois = s a P L bp ■ spP L s a , (3) 

O23 = SaCr^Pjiba ■ spa^PRSp, O 2 4, = s a a fn 'P R b l 3-S0a lll ,P R s a , (4) 
O25 = s a a^P L b a ■ spa^P L sp, 2 6 = s Q cr A " y P L bp ■ spa^P L s a , (5) 

satisfies the polarization Anomaly Resolution Criteria (pARC) [9], allowing one to possibly solve the polarization 
anomaly in B°(B°) -> (j)K*°(K*°)-<lec&y. Here, a and J3 are the color indices. In what follows we denote the set of 
operators in Eqs. ©-© as pARC operators. In Ref. 0] it was noted that the (pseudo-)scalar operators O15-18 can 
be expressed in the basis of (pseudo-)tensor operators O23-26 by Fierz transformation 

Oi 8 = ^024 -^23, (6) 
Ois = ^0 26 - \o 25 . (7) 
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The contributions of the operators Ois_26 to the helicity amplitudes of B° — > (j)K*°-dec&y can be calculated within 
the QCD factorization (QCDF) approach in terms of the corresponding Wilson coefficients and hadronic form factors. 
In this approach the helicity amplitudes take the form 0, [§] 



H 00 = -4ifJm 2 B (d 23 - d 25 ) \h 2 T 2 (m\) - ft, 3 T 3 (m|)] , 

H±± = -Ufjml {a 23 [±fxT x {ml) - / 2 T 2 (m|)] + a 25 [i/iT^m 2 ^ + / a T 2 (mJ)] } . 
Here the </>-meson tensor decay constant fj and the form factors of the B — K* transition are defined as 

{K*{p\e 2 )\sa^{l+ lb ) S \B{p)) = 2te^ pt7 e 2 '*prp'°T 1 (q 2 ) + T 2 (q 2 )[e^(m% - m\,) 



(e* 2 - P )(p + p'),}+T 3 (q 2 )(e* 2 -p) 



q 2 {p + p')v 



K* 



(8) 
(9) 

(10) 

(11) 



with q = p — p' and m B = 5.279 GeV, m K * — 0.892 GeV, — 1.019 GeV being the masses of the B 7 K* and 
mesons, respectively. The kinematical factors in Eqs. ((H)) are: 
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(12) 
(13) 
(14) 



where p c is the momentum of the <j> or K* meson in the rest frame of the decaying B° meson. The effective coefficients 
in Eq. ([5]) are expressed in terms of the Wilson coefficients as @ 



a 23 = yl 
a 25 = 1 
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(15) 
(16) 



The last terms in ([T5|) and p^|) indicate corrections due to deviations from the QCDF. 

On the basis of the above equations in Ref. [9( there has been made an analysis of the experimental data obtained 
by BABAR 2] and BELLE |f on the angular distribution in B°(B°) -> c/)K* Q (K*°)-dec&y. It was shown that there 
are two theoretical scenarios which can separately account for the polarization anomaly of these data. 
Scenario (i): a 23 = and 



|a 36 | = 2.10lg;Jjj x 10" 4 , S 25 = 1.15 ± 0.09, 



Scenario (ii): 025 = and 



IS23I = i.7oig;S x 10" 



S 23 = 2.36 ±0.10, 



-0.12 ±0.09. 



= 0.14 ±0.09. 



(17) 



(18) 



Here the following notations were used dij = \dij\e ij e 1 ^ , identifying 0y and <Jy with the weak (coming from the 
terms in Eq. |T])) and strong phases, respectively. The values of Eqs. (fTTj) and ([15]) correspond to the best fit values 
for the combined data of BABAR and BELLE [|. In what follows we use these results as a criterion to assess if 
a particular model is able to resolve the polarization anomaly in question or not. 



III. PARC OPERATORS IN fl p SUSY 

In relation to the polarization anomaly in B° — > (j>K*° we are studying the AB = 1 transitions on the quark level. 
Here we derive the effective Lagrangian describing these transitions within the minimal $ p SUSY model ($ p MSSM) 
and show that among the resulting set of operators there emerge the pARC operators O15 and On. In the generic 
case of Ifap MSSM R-parity is violated by the following terms in the superpotential 

W $p = fijLjH 2 + ^uUL.El + K^QM + ^X'l^UfD^Dl, (19) 
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and in the soft SUSY breaking part of the scalar potential 

Vg lt = ^kUL 3 El + K.L.Qjr, + PJ^UtDpt + J% 3 ~L d H 2 + ~^L 3 H\ + H.c. (20) 



In Eq. (fl9|) L, Q stand for the lepton and quark doublet left-handed superfields, while E c 1 U c , D c for the lepton and 
itp, down quark singlet superfields; H.% is the Higgs doublet superfields with a weak hyperchargc Y — +1, respectively. 
In Eq. (|20|) Li denotes the scalar sleptqn weak doublet, H1.2 are the scalar Higgs doublet fields. In the above equations 
the trilinear terms proportional to A, A', A, A' and the bilinear terms violate lepton number, while the trilinear terms 
proportional to A", A" violate baryon number conservation. The coupling constants A (A") are antisymmetric in the 
first (last) two indices. The bar sign in A', A" denotes that all the definitions are given in the gauge basis for the quark 
fields. Later on we will change to the mass basis and drop the bars. Using the freedom in the definition of lepton and 
Higgs superfields we choose the basis where the vacuum expectation values of all the sneutrino fields vanish: (vi) — 0. 

The Lagrangian terms generated by the trilinear terms of the superpotential in Eq. (|19p and involving two down 
quarks needed for the construction of the pARC operators in (H])-© are as follows: 

£a = -KjkViLdkPhdj - ^"jkUi R djPLd c k + H.c. (21) 

where dj stands for the down quark. It can be easily seen that the interactions in Eq. (|21[) can generate in second 
order perturbation theory the only AB = 1 effective operator contributing to B°(B°) — » 4>K*° (K*°)-deca.y. This is 
the operator (sPLb)(sPus) which does not belong to the pARC operators listed in Eqs. (H])-©. Thus, we conclude 
that the trilinear ^-couplings alone cannot resolve the polarization anomaly in B°(B°) — > <fiK*° (K*°)-deca.y. 

Let us see if the bilinear -terms may help in the solution of this problem. The presence of the bilinear terms 
leads to terms in the scalar potential which are linear in the sneutrino fields, Vi. First, this results in v — H® 2 
mixing. Also, the linear terms drive the Vi fields to non-zero vacuum expectation values (vi) 7^ at the minimum 
of the scalar potential. At this ground state the MSSM vertices Zv v and Wev produce the gaugino-lepton mixing 
mass terms Zv(v) : We(v) (with W, Z being wino and zino fields). These terms taken along with the lepton-higgsino 
HiLiHi mixing from the superpotential of Eq. (|19[) form 7x7 neutral fermion and 5x5 charged fermion mass 
matrices (ll| . This leads to a non-trivial neutrino mass matrix and Lepton Flavor Violation in the sector of charged 
leptons. However, these effects are obviously irrelevant for the generation of the effective 4-quark operators. 

The above mentioned effect of sneutrino-Higgs mixing v — H\ 2 is different. It corresponds to a non-diagonal mass 
matrix for the neutral scalars (H\ , H® , v ei i>^,v T ) in the bilinear part of the $ p scalar potential [l2j]. From Eqs. (jTT 
and (12011 we write: 



y% = (^* H H\ + fa H 2 + filjH^Lj + H.c. (22) 

Using the minimization condition 

filj + fi*fij + tan/3 = (23) 
in the basis of lepton and Higgs superfields where (vi) = we can rewrite Eq. (|2"2")l in the form 



^° ft = Aaj (#2 - tan PH{)Lj + H.c. , (24) 

where tan/3 = (-fff) / (Hi)- Rotating these fields to the mass eigenstate basis we assume smallness of sneutrino-Higgs 
mixing characterized by the small ratio (jlkj /Mj ll 2 ) 2 , where /Ijy is the ^ p soft parameter from Eq. (|20[) and 2 
are the neutral Higgs masses [l3j . In the leading order in this small parameter we obtain the following interactions of 
sneutrinos with down quarks and charged leptons 



Cm = rjj 

with the couplings 



-T7—{<k di) + -jrj—(li U 



(25) 



Vj = ^7rrm[^l~^ ] - (26) 



gi _ 2 tan/3 / cos a 

Here a is the mixing angle of the neutral Higgses in the limit of no mixing with the sneutrino fields: 

B.\ =-sina-/i? + coso;-/i2, = cosa • h\ + sina • h% , (27) 
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where h® 2 are the corresponding mass eigenstates with the masses Mh t , Mh 2 - Note that , which has no couplings 
to the down quarks and leptons, does not contribute to Eq. (|2"5l) . 

Now, combining the trilinear and bilinear ]fi p -interactions from Eq. (|2Tj) and Eq. (|25|) . as shown in Fig.l, we obtain 
in second order perturbation theory the following effective Hamiltonian after integrating out the heavy sneutrino 
fields: 

= ^r(d 3 d 3 ) (-^K: n3 d m p R b + ^K 3m d m p L b 

1V1 W \" l Ui " l i>i 

+ ^rih h) (-\K%a d m P R b + -J-A< 3m d m P L b) + H.c. (28) 

The 4-quark terms involve the pARC operators O15 and On from the list of Eqs. ([I])-© with the following Wilson 
coefficients: 



Thus ]jtp SUSY seems to satisfy the pARC as it allows appropriate operator structures. In the following we have 
to check if the existing experimental constraints on the -parameters entering into the definition of the Wilson 
coefficients allow one to accommodate the values of Eqs. (TT7|) and (fT5|) . 



IV. EXPERIMENTAL CONSTRAINTS ON WILSON COEFFICIENTS 



Examining Eq. (|28[) we note that the strength of both the 4-quark and quark-lcpton operators is determined by the 
same combination of the R-parity conserving and ^-parameters forming the Wilson coefficients 015,17. Therefore, one 
can directly constrain the C15 17 parameters from the existing stringent experimental upper bound on the B s 
branching ratio [l5T | 



Br(B s -> < 1.0 x 1(T 7 (90% C.L.) . 



(30) 



An important advantage of this constraint is that it applies to the coefficients 015,17 as a whole, avoiding uncertainties 
related to the presence of several R-parity conserving (tan/3, a, M^ ^jinp) and violating parameters (/xL, A'). 



The contribution of the quark-lepton interactions in the Lagrangian 
written in terms of the Wilson coefficients C15 17 as 



to the decay rate of this process can be 



r(B s 



2 32tt 



m s J \ mb 



D , 



(Cl 5 - Ci 7 ) 2 



where we used 



2m p 



3/2 



(31) 



\s l5 b\B° s ) = ih 



mt, + m s 



(32) 



We use the following numerical values for the quantities in above equations: Jb s — 0.2 GeV [l6|, rriB s = 5.367 GeV, 
nib = 4.6 GeV, m s — 0.15 GeV and tb 3 = 1-46 x 10~ 12 s [l|. Considering the two scenarios of Ref. as displayed in 
Eqs. (fTT]) and (fT8|) [denoted as (i) and (ii)] we get from the experimental limit (|30|) the following upper bounds 

|ci 5 |,|ci 7 | < 1.4 x 10- 4 . (33) 

Using the definitions of Eqs. (fTS"|) and (fTB|) these limits can be translated to upper limits on the effective coefficients 



laaaUassI < 5.9 x 10" 



(34) 



These limits are about two orders of magnitude smaller then the values given in Eqs. (fTT)) and (jT5J) required for the 
solution of the polarization anomaly. 

Thus, we conclude that the polarization anomaly observed in B°(B°) — > <j)K*°(K*°) decay by the BABAR Q and 
BELLE 3 collaborations cannot be explained within the I/tp SUSY framework, despite the occurrence of effective 
operators with the chiral structure required qualitatively. 
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As a byproduct of our analysis the limits of Eq. ([33|) set new upper limits on the products of the soft and super- 
potential Ifip -parameters of Eqs. (p?l)|) and (|29[) . Since the expressions for the Wilson coefficients 015,17 contain the 
R-parity conserving parameters as well we choose one representative point in the SUSY parameter space in order to 
illustrate the limits on the -parameters. We take a typical mSUGRA: the so-called SPS la point from the list of 
nine Snowmass benchmark points [l7j]. This choice corresponds to tan/3 = 10, mo = — = Q.2hm l / 2 = 100 GeV 
and (x > 0. For this parameters we find 

To our knowledge in the literature (for a review see, for instance [14|) there have not been established experimental 
limits on these products of ljt p -parameters. However, there exist bounds on fi^, AJ23 and X' i32 separately from various 
low energy processes flij . This allows one to obtain indirect bounds on their products and compare them with 
those in Eq. (|35|) . The soft $ p -parameter /Lt 2 ^., contributes to the neutrino mass matrix at one-loop level. Thus it is 
constrained by the present limits on neutrino masses and mixing from neutrino oscillations. With the SPS la set of the 
R-parity conserving parameters one has: (/i2i/100 GeV) 2 < 10 -4 . Existing constraints on the trilinear ^-couplings 
are typically as follows: A^ 2 3, X' i32 < 0.2. Combining these constraints we have the limits: 

(i^) 3 |^l'(l5ifev) 2 W-l £2 - 0><10 " <36 » 

which are two orders of magnitude better then those in Eq. Nevertheless, the latter can still be useful as direct 

constraints on the specific products of the bilinear and trilinear ^.-parameters. Note that these constraints correspond 
to a particular point in the MSSM parameter space and in some other points the above limits may significantly change. 
The detailed study of this question is beyond the scope of the present paper. 

V. CONCLUSIONS 

We analyzed the I/l p SUSY model with respect to its ability to account for the polarization anomaly in B°(B°) — > 
cj)K* (K* )-decay observed by the BABAR Q and BELLE Q collaborations. Within this framework we have deter- 
mined the effective AB — 1 operators with chirality structures appropriate for a possible resolution of this anomaly. 
However, the experimental data on B — ► /i + /i - -decay set stringent limits on the respective Wilson coefficients, which 
are about two orders of magnitude below the values required to resolve the polarization anomaly. This gap of two 
orders of magnitude can hardly by eliminated by the uncertainties in the hadronic parameters involved in the calcu- 
lation of the helicity amplitudes of B°(B°) — > <^?r* (.K"* )-decay. Therefore, we do not believe that I/L p SUSY is able 
to account for the B-decay polarization anomaly. 

As a byproduct we used the experimental data on B — > fi + \i~ -decay to set a new upper limit on the product 
of the two Jfi p -parameters Jj^ | A^ 23 1 and pfe l-Maa I) where and X' i j k are bilinear soft and trilinear superpotential 
I/tp -parameters, respectively. 
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FIG. 1: The Ijt p SUSY contribution to the b — > sss (a) and to the bs —> 11 transition operators. The sign ® denotes $ p soft 
sneutrino-Higgs mixing. The left hand vertices in both diagrams are due to the fy v superpotential A' coupling, while the right 
hand ones correspond to the R-parity conserving Hi — q — q and Hi — I — T Yukawa couplings. 
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